Quantification of noise in the bi-functionality induced post-translational modification 
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We present a generic analytical scheme for the quantification of fluctuations due to bi-functionality 
induced signal transduction within the members of bacterial two-component system. The proposed 
model takes into account post-translational modifications in terms of elementary phosphotransfer 
kinetics. Sources of fluctuations due to autophosphorylation, kinase and phosphatase activity of 
the sensor kinase have been considered in the model via Langevin equations, which are then solved 
exactly within the framework of linear noise approximation. The resultant analytical expression 
of phosphorylated response regulators are then used to quantify the noise profile of biologically 
motivated single and branched pathways. Enhancement and reduction of noise in terms of extra 
phosphate outflux and influx, respectively, have been analyzed for the branched system. 

PACS numbers: 87.18.Mp, 87.18.Tt, 87.18.Vf 



I. INTRODUCTION 

The response of living systems in presence of an exter- 
nal stimulus is coordinated by highly specialized signal 
transduction machinery. In the bacterial kingdom this 
is achieved by the well characterized two-component sys- 
tem (TCS) minimally comprised of the membrane bound 
sensor kinase (SK) and the cytoplasmic response regula- 
tor (RR) [1-3]. The machinery of TCS is utilized by 
the bacteria to process the information of external sig- 
nals in terms of phosphotransfer kinetics. When applied, 
an external stimulus causes phosphorylation at the his- 
tidine residue of SK which then gets transferred to the 
cognate (and/or non-cognate) RR at their aspartate do- 
main. The phosphorylated RR then acts as a transcrip- 
tion factor for several downstream genes, as well as for 
the activation/represion of its own operon. It is now a 
well established fact that in addition to being a source 
(kinase), a SK can equally act as a sink (phosphatase) 
while interacting with a RR. Such bi-functional behavior 
of a SK towards RR can altogether build a robust motif 
in the bacterial signal transduction network [ij-Q . 

Expression of proteins in individual cell is usually 
driven by the fluctuations present within the cellular en- 
vironment as well as the fluctuations imposed by the ex- 
ternal stimulus This often leads to variability in 
the expression level within the context of single cell [1 5l — 
[l8j . When observed in the bulk such fluctuations gets 
averaged out over the cellular population. The prevalent 
fluctuations, whether external or internal, do not only 
effect the dynamics of gene expression it pla ys a major 
role in post-translational modification as well |19l . |20| . In 
this connection it is also important to mention the role 
of cellular fluctuations in the different signal transduc- 
tion motif that primarily uses phosphotransfer mecha- 



nism. Using a push-pull amplifier loop mechanism theo- 
retical study has been made to analyze the signal trans- 
duction within the photoreceptor of retina j2l[ . Theoret- 
ical model has been proposed to study the effect of re- 
versibility in the phosphorylation-dephosphorylation cy- 
cle that can generate bistable behavior in the presence 
of noise and can propagate within the signaling cascade 
[22I ]. In the context of robustness in bacterial chemotaxis, 
reversibility on a signaling cascade has been shown to ex- 
ert a stabilizing effect of adaptation through methylation 
[23| . Correlation between extrinsic and intrinsic noise 
that arises due to external signal and internal biochemi- 
cal pathways, respectively, has been reported to enhance 
the robustness of zero-order ultrasensitivity (2J| . 

Post-translational modification in terms of phosphate 
transfer is important to generate the pool of phospho- 
rylated RR that act as a transcription factor for several 
downstream genes. Bi-functionality on the other hand 
plays a crucial role to maintain this pool as the informa- 
tion flows through the phosphotransfer motif. Thus bi- 
functionality and post-translational modification works 
hand in hand to maintain the optimal pool of phos- 
phorylated RR. Since this composite functional behavior 
take place in a noisy cellular environment it is worth- 
while to investigate the role of cellular noise on the bi- 
functionality controlled post-translational modification 
of the components of the well composed signal transduc- 
tion machinery. The specific question we aim to address 
here is how the cellular fluctuations affect this functional 
motif or in a noisy environment whether individual cell 
can maintain bi-functionality while generating optimal 
level of noise? To address this questions we develop 
in the present communication a generic mathematical 
model to study the role of molecular noise on the post- 
translational modification of system components and on 
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FIG. 1. (color online) The wiring diagrams for the proposed 
post-translational interactions between SK and RR. The cyan 
and magenta blob stands for the SK and RR, respectively, (a) 
The generalized m:n system with cognate (solid black arrow) 
and non-cognate (dashed colored arrow) kinase and phos- 
phatase interactions, (b), (c) and (d) are wiring diagram for 
1:1, 1:2 and 2:1 interactions, respectively. The boxed wiring 
diagram depicts the phosphate addition (+P) and removal (- 
P) kinetics between a pair of cognate and/or non-cognate SK 
and RR that ultimately results into phosphorylated RR (ma- 
genta blob with a yellow blob on top). The dotted arrowhead 
is for auto-dephosphorylation of RR. For simplicity we do not 
show the yellow blob in (a)-(d). 



the bi-functionality of the signal transduction machinery 
of bacterial TCS. 



II. THE MODEL 

We start by considering a simple system describing 
post-translational modification due to phosphotransfer 
mechanism of a typical TCS, where m numbers of SK 
interacts with n numbers of RR, the ultimate product of 
which is R p , the phosphorylated RR. We call the pro- 
posed model as the mm system (Fig. QJi) where each of 
the SK, RR and their phosphorylated forms are desig- 
nated as S, R, S p and R p , respectively. The generic 
model we consider here involves single pair interaction 
(Fig. [Tp). In addition it takes care of branched pathways 
[3|; for example the 1:2 system (Fig. [TJ;) mimics the one 
to many pathway as observed in chemotaxis system in E. 
coli, where the SK CheA phosphorylates two RR, CheY 
and CheB (25[. Similarly the 2:1 system (Fig. [Hi) follows 
the kinetics of many to one pathway as observed in V. 
cholerae where the SK LuxS and CqsS phosphorylate the 
RR LuxO (26|, respectively. 

As mentioned earlier, in typical bacterial TCS the key 
steps of phosphotransfer mechanism involve autophos- 
phorylation at SK, transfer of phosphate group from SK 
to RR, and SK mediated removal of phosphate group 
from RR (see boxed diagram in Fig. Q}. To keep the 
model simple we do not consider the synthesis (birth) 
or degradation (death) of any system component. The 
interaction we consider here may be of cognate and/or 



non-cognate type. Out of the mm pair one can consider 
the specific interaction between i-th SK and j-th RR, 
where 1 ^ i ^ m and 1 ^ j $J n, to write down the 
elementary kinetic steps considering the minimal inter- 
action between a specific pair 
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In the above kinetic steps, Eq. (Hal) considers autophos- 
phorylation at the histidine residue of the SK. Generally, 
autophosphorylation takes place under the influence of 
an external signal [J-Q which we consider to be of con- 
stant type and have been absorbed in the rate constant 
ki. Eqs. (II btilc[) takes into account the kinase and phos- 
phatase activity of the SK, respectively, thus considering 
the bi-functional behavior of the SK. Eq. (jldl) denotes 
the auto-dephosphorylation of the RR independent of the 
phosphatase effect of SK on RR [6| . 

Due to the inherent noisy nature of the cellular en- 
vironment each of the four reactions mentioned above 
are influenced by fluctuations and in turn affect the copy 
numbers of each system component. To take this into 
account we introduce Langevin noise terms that can 
influence each of the reactions independently given by 
Eqs. (|lallld[) . The interaction of a single SK with multiple 
RR or vice versa in presence of fluctuations we consider 
here can be compared with stochastic system-reservoir 
formalism where a single system interacts with multi- 



ple reservoir or vice versa [27] . In terms of phosphate 
transfer, SK pumps in energy to RR when act as kinase, 
whereas the same withdraws energy from RR while act- 
ing as phosphatase. The stochastic differential equations 
describing the phosphorylated SK and RR in presence of 
fluctuations can be written as 

dS ■ n 
^ — ki(Sxi ~ Spi) — k_iSpi — ^ ^ kkijSpi 

3=1 

x(R Tj - Rpj) + &(t) (2a) 



dRp. 
dt 



fcfeij Spi (Rtj Rpj ) 2_j ^ pi 3 

i=l i=l 

x (S Ti - Sn) + Kp^Rpj + tpj (t). (2b) 



Here Sri — Si + S P i and Rxj — Rj + Rpj stand for the 
total amount of i-th SK and j'-th RR, respectively. The 
noise terms £j and ipj take care of the fluctuations in 
the copy number of Spi and Rpj, respectively. Within 
the framework of linear noise approximation we define 
the statistical properties of the Langevin terms obeying 
the fluctuation-dissipation relation 0,128-30] with zero 
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mean, (£t(i)) = (ipj(t)) — and 

(&(*)&(* + t)> = 2fc 2 (fe - 5 P4 )<5(r), 

m 

(ipjiWjit + t)> = 2^ kkijSpiiRrj - Rpj)S(j), 

i=l 

with (Spj) and (R P j) being the mean values at the steady 
state. In addition, we consider both the noise terms are 
correlated [24, .311 



(ZiitfTpjit + T )) = -k kij S pi {R Tj - R pj )5(r). 

Linearizing Eqs. (|2all2bp around the steady state, i.e., 
S p i = (S p i) + SSpi and R p j = (R P j) + SR P j, we have 
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Solving Eq. ([3]) and performing Fourier transformation 
= SX (t)e~~ lut dt on the resultant solution we 
have in matrix notation the generalized solution for both 
6S p (ui) and 5R p (uj), 



6S p (lu) = A- 1 (<S pK )<5R P M + |( w ; 
ffip(w) = P 1 (b t A- 1 £(w) + ■$(<*]) 



(4a) 
(4b) 



where P = C — B t A _1 (S p k)- In the above expressions 
(|4all4b[) , dS p and STl p are mxl and nxl dimensional col- 
umn vectors with elements SS P i and SR P j, respectively. 
Likewise, £ and if) are mxl and nxl column vectors 
with elements and t/>j-, respectively. A and C are di- 
agonal matrix of order m x m and n x n with elements 
(iw + aj) and (iw + Cj), respectively. Additionally, (S p k) 
and B are matrix of order m x n with elements (S P i)kkij 
and 6^ , respectively. 

Since we are interested in the effect of noise on phos- 
phorylated RR, R p , which acts as transcription factor for 
one or more genes including the gene that encodes SK 
and RR we now focus on the solution of Eq. (|4b|) only. 
From the structure of Eq. (|4bl) it is clear that dynamics of 
Rp is now decoupled from S p . To understand the role of 
fluctuations in phosphotransfer processes, we define the 
quantity noise at steady state, t]ji p = gr /(Rp), While 
calculating the noise for the three different systems (1:1, 
1:2 and 2:1) mentioned in Fig. (flj we only focus on the 
noise level of R p \ , the phosphorylated form of R\ . Noise 
generated due to other interactions (Si and i?2, and S2 
and Ri ) are considered to add extra layers of information 
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FIG. 2. Plot of noise, steady state protein level, time se- 
ries and power spectra for 1:1 system, (a) Noise r/R pl as a 
function of log 10 kii. The noise profile has been shown to con- 
tain three different regions, viz mono-functional region (fen < 
feapi)) crossover region (fcn w k ap i) and bi-functional region 
(fcn > k ap i). In the absence of auto-dephosphorylation kinet- 
ics (Eq. (|ldl) ~) only the bi-functional domain exists, whereas in 
the absence of phosphatase kinetics (Eq. (JTcJ ) only the mono- 
functional domain becomes prevalent, (b) Steady state R p i 
as a function of log 10 feu. (c)-(e) Time series of R p i for low 
(fcn = 1CP 5 ), intermediate (fen = 10 -2 ) and high (fen = 1) 
values of feu, respectively, generated using Gillespie algorithm 
H3;H3| (f) Normalized power spectra for low (solid line), in- 
termediate (dashed line) and high (dotted line) values of fen. 
In all the cases fci/fc_i = 5, k ap i = 0.01 and Sti ~ Rti ~ 20. 



on the noise profile of R p i . During the calculation of noise 
and power spectra we have considered kkij — k p ij = kij 
(i,j = 1,2), as kinase and phosphatase activity between 
SK and RR takes place via protein-protein interaction. It 
is important to note that while interacting with its part- 
ner a SK shows both mono-functional and bi-functional 
behavior for k ap j > ka and k a ni < 

At knni ~ ki 



^apj ^ <n 3 ^ apj ^ n^j, respectively [6|. 
^ apj ~ mj (cross over regime) system makes transition 
from mono- to bi-functional domain. 



III. RESULTS AND DISCUSSION 

In Fig. (fJk) noise profile of R p i has been shown in a 
semilog plot. For 1:1 system at a low value of fen noise 
has a nonzero value which goes down as fen value in- 
creases. As fen value increases further noise increases and 
reaches a higher value. As evident from the definition, 
noise is inversely proportional to the population of steady 
state Rpi. To check the role of (Rpi) on steady state 
noise we have calculated (Rpi) as a function of log 10 fen 
(Fig. fJb) from which it is clear that the protein profile 
develops exactly in the opposite way of the noise profile 
and imparts a inverse effect on the noise development. 
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FIG. 3. Plot of noise and steady state protein level for 
branched (1:2 and 2:1) pathways. (a,c) Noise profile for 1:2 
and 2:1 system, respectively, as a function of log 10 fcn. The 
solid, dashed and dotted lines are for low (10 ), intermediate 
(10 -2 ) and high (1) values of k 12 (for 1:2) and k 2 i (for 2:1). 
(b,d) Steady state level of R p i for the same values of fci2 and 
&2i as in (a,c). For all the cases ki/k-i = 5, k ap j — 0.01 and 
Sn = Rt 3 = 20. 



For the low value of kn the auto-dephosphorylation k ap \ 
dominates over the phosphatase activity of SK on RR 
(fen < k ap i). In this regime the sensor shows mono- 
functional activity by acting as a kinase only, which is 
still lower than k ap \. This effectively reduces the R p \ 
level (Fig. and increases the noise of the system. In 
the limit kn ~ k ap i (vertical dotted line in Fig. [5^) R p i 
level reaches at its maximum (Fig. [2J5 and Fig. [2]i) while 
reducing the noise. When kn exceeds k ap i (kn > k ap i), 
the phosphatase activity of SK starts to show up in ad- 
dition to its kinase activity. In this regime bi-functional 
property of SK comes into play reducing the copy of R p i 
(Fig. henceforth increasing the noise of the system. 
To understand how the system relaxes under the influ- 
ence of the protein protein interaction, we calculate the 
power spectra S(ui) = (SR p i(uj)SR p1 (u}')) (Fig. Off). The 
resultant spectral lines are plotted for low, intermediate 
and high kn (= kkii = k p n) values. As expected the 
power spectra relaxes faster for low kn value compared 
to intermediate kn value which again relaxes faster com- 
pared to higher kn value. For a low value of kn the 
conversion of R into R p is a slow process and hence fast 
fluctuations in the copy number have minimal effect on 
the power spectrum. As kn value increases the conver- 
sion rate increases and thus gets affected by the noise in 
the copy number which results into slower relaxation. 

In Figs. (|3k-[2t) hi a semilog plot we show the noise, 
rin 1 for 1:2 and 2:1 system as a function of kn- For 
comparison we refer to noise profile of 1:1 system shown 
in Fig. (Hk). In 1:2 system a single SK, Si interacts 



with two RR, Ri and R2, with its bifunctional properties 
acting on both the RRs. In Fig. (J5Ji) we show the noise 
generated for R p i while taking into account the kinase 
and phosphatase rates (kki2 — fcpi2 = &12) between Si 
and i?2 to be low (ki2 = 1CP 5 ), intermediate (fci2 = 
1CP 2 ) and high (ki2 — 1). For low and intermediate 
fei2 value the noise profile looks almost like 1:1 system 
as fcn is varied. This happens as interaction between 
Si and i?2 adds a weak layer of information on Ri due 
to mono- functional property of Si on R2 (k ap 2 ^12)- 
Whereas for high ki2 value a huge amplification of noise 
occurs (the dotted line in Fig. In this domain as 

k ap 2 < ki2, SK starts to show its bi-functional property 
and is more active in its interaction with R2 , rather than 
with i?i . Such active interaction between Si and R2 adds 
an extra layer of outflux of phosphate group from Ri (the 
dotted line in Fig. [3Jd) thus leading to low level of (Rpi) 
and enhancement of noise. 

In 2:1 system a single RR, i?i interacts with two SK, 
Si and 5 2 - In Fig. (J3J:) we show the noise generated for 
R p i while taking into account the kinase and phosphatase 
rates (fcfc 2 i = k P 2i = &2i) between 5 2 and Ri to be low 
(k 2 i = 10~ 5 ), intermediate (k 2 i = 10~ 2 ) and high (fc 2 i = 
1). For low value of A:2i the noise profile is similar to 
that one for 1:1 system as kn is increased. Although 
in this domain S2 acts as kinase only, it provides a low 
level of input on Ri as k ap i > /c2i- Interesting behavior 
emerges as £21 takes intermediate and high value. In the 
intermediate domain maximal level of (Rpi) is produced 
due to extra influx of phosphate group. This large influx 
due to fc 2 i can overpower the low kinase effect of kn, 
henceforth increases the steady state level of R p i as a 
effect of which the noise level attains a minimum value. 
For high fc 2 ij k ap i < fei where S2 starts to show its 
bi-functional property via phosphate input and removal. 
This helps the composite system to maintain a high value 
of noise for a wide range. Note that, compared to the low 
and intermediate domain the protein level in this region 
goes down drastically due to strong phosphatase activity 
of S2 on i?i . It is interesting to note that for intermediate 
and high k2i value the composite system looses it mono- 
functional behavior almost completely (Fig. [5Ji) . 



IV. CONCLUSION 

To conclude, we have provided a generic description for 
the calculation of noise due to post-translational modi- 
fication in bacterial TCS. From exact analytical calcula- 
tion within the purview of linear noise approximation it 
is possible to quantify the steady state noise for a sin- 
gle pair and branched pathways. For single pair system 
our analysis shows the effect of bi-functionality of SK on 
noise generation and could differentiate the mono- and 
bi-functional domain in the noise profile. Calculation for 
branched pathways shows enhancement and reduction of 
noise for the composite system in terms of extra phos- 
phate outflux and influx, respectively. Our analysis sug- 
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gests that in one to many system as in the chemotaxis 
pathway of E. coli enhancement of fluctuations happens 
due to extra outfiux of phosphate group within the mem- 
bers of TCS. On the other case, for many to one system 
mimicking the quorum sensing network of V. cholerae an 
optimal level of noise can be maintained via extra influx 
of phosphate group. To maintain such low noise activity 
SK of V. cholerae phosphotransfer circuit might prefer 
to operate in the cross over domain. 
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